GEOMETRY AND COMPLEXITY OF O'HARA'S ALGORITHM 



MATJAZ KONVALINKA AND IGOR PAK 

Abstract. In this paper we analyze O'Hara's partition bijection. We present three type 
of results. First, we show that O'Hara's bijection can be viewed geometrically as a certain 
scissor congruence type result. Second, we obtain a number of new complexity bounds, 
proving that O'Hara's bijection is efficient in several special cases and mildly exponential 
in general. Finally, we prove that for identities with finite support, the map of the O'Hara's 
bijection can be computed in polynomial time, i.e. much more efficiently than by O'Hara's 
construction. 



1. Introduction 

Ever since the pioneer work by Sylvester and his school, there has been a quest to find 
bijective proofs of many interesting partition identities. Despite remarkable successes in 
the last century (see |P3j ) and some recent work of both positive and negative nature (see 
e.g. \P2\ IP4| ). the problem remains ambiguous and largely unresolved. Much of this stems 
from the lack of clarity as to what exactly constitutes a bijective proof. Depending on 
whether one accentuates simplicity, ability to generalize, the time complexity, geometric 
structure, or asymptotic stability, different answers tend to emerge. 

In one direction, the subject of partition bijections was revolutionized by Garsia and Milne 
with their involution principle [GMll IGM2] . This is a combinatorial construction which 
allows to use a few basic bijections and involutions to build more involved combinatorial 
maps. As a consequence, one can start with a reasonable analytic proof of a partition 
identity and trace every step to obtain a (possibly extremely complicated) bijective con- 
struction. Garsia and Milne used this route to obtain a long sought bijection proving the 
Rogers- Ramanuj an identities, resolving an old problem in this sense [GM2j . Unfortunately, 
this bijection is too complex to be analyzed and has yet to lead to new Rogers- Ramanuj an 
type partition identities. 

After Garsia-Milne paper, there has been a flurry of activity to obtain synthetic bijections 
for large classes of partition identities. Most of these bijections did not seem to lead any- 
where with one notable exception. Remmel and Gordon found (rather involved) bijective 
proofs of basic partition identities due to Andrews [Rl [G] . The latter are direct extensions 
of Euler's distinct/odd theorem and have a similar straightforward analytic proof [X t IP3| . 
Then O'Hara made a surprising discovery that Remmel's and Gordon's bijections can be 
streamlined to give the same bijective map with a simple construction |Q1[ I02| . In fact, 
O'Hara proved that the resulting bijection is a direct generalization of Glaisher's classi- 
cal bijection proving Euler's theorem. Moreover, in her thesis [01], O'Hara showed that 
her bijection is computationally efficient in certain special cases. Until now, the reason 



Date: February 2, 2008. 



2 



MATJAZ KONVALINKA AND IGOR PAK 



why O'Hara's bijection has a number of nice properties distinguishing it from the other 
"involution principle bijections" remained mysterious. 

In this paper we obtain results of both positive and negative type. First, we analyze the 
complexity of O'Hara's bijection, which we view as a discrete algorithm. We prove a general 
result (Theorem l3.2p . which given an exact formula for the number of steps of the algorithm 
in certain cases. From here we show that O'Hara's bijection is computationally efficient in 
many special cases. On the other hand, perhaps surprisingly, we prove that the number of 
steps can be (mildly) exponential in the worst case (Theorem 13.71 part In fact, even 

when the natural speed-up is applied, the worst complexity does not improve significantly 
(see Section [5]). This is the first negative result of this kind, proving the analogue of a 
conjecture that remains open for the Garsia-Milne's "Rogers- Ramanuj an bijection" (see 
Subsection 16. ip . 

Second, we show that O'Hara's bijection has a rich underlying geometry. In a manner 
similar to that in [PTl IPV| . we view this bijection as a map between integer points in 
polytopes which preserves certain linear functionals. We present an advanced generalization 
of Andrews's result and of O'Hara's bijection in this geometric setting. In a special case, 
the working of the map corresponds to the Euclid algorithm and, more generally, to terms 
in the continuing fractions. Thus one can also think of our generalization as a version of 
multidimensional continuing fractions. 

Finally, by combining the geometric and complexity ideas we show that in the finite di- 
mensional case the map defined by O'Hara's bijection is a solution of an integer linear 
programming problem. This implies that the map defined by the bijection can be com- 
puted in polynomial time, i.e. much more efficiently than by O'Hara's bijection. This 
suggests that perhaps in the general case O'Hara's bijection can also be sped up to work in 
polynomial time (see Subsection 13. Sp . 

The paper is structured as follows. We start with definitions and notations in Section [2l In 
Section [3l we describe the main results on both geometry and complexity. Proofs of most 
results are postponed until Section [H We conclude with a quick application in Section [5] 
and final remarks in Section [H 



2. Definitions and background 

2.1. Andrews's theorem. A partition A is an integer sequence (Ai, A2, . . . , A^) such that 
Ai > A2 > . . . > A^ > 0, where the integers Aj are called the parts of the partition. The sum 
n = Aj is called the size of A, denoted |A|; in this case we say that A is a partition of 

n, and write A h n. We can also write A = l"*i2'"2 . . . ^ where rrii = mi{\) is the number of 
parts of A equal to i. The support of A = l™i2'"2 ... is the set {i: rrii > 0}. The set of all 
positive integers will be denoted by P. 

Denote the set of all partitions by V and the set of all partitions of n by Vn- The number 
of partitions of n is given by Euler's formula 

00 00 ^ 
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For a sequence a = (ai, a2, . . .) with Oj G P U {oo}, define A to be the set of partitions A 
with mi{X) < ai for all i; write An = ACiVn- Denote by supp(a) = {i: < 00} the support 
of the sequence a. 

Let a = (oi, a2, . . .) and b = 62, . . .)• We say that a and b are ip- equivalent, a b, if (p 
is a bijection supp(o) supp(&) such that iai = (p{i)b^(^i-) for all i. If a 6 for some if, we 
say that a and 6 are equivalent, and write a ^ b. 

Theorem 2.1 (Andrews) Ifa^b, then \ An\ = \Bn\ for all n. 

Proof. We use the notation = 0. Clearly, 

n=0 i=l j=l n=0 

which means that \An\ = □ 

Consider the classical Euler's theorem on partitions into distinct and odd parts. For a = 
(2, 2, . . .) and b = (00, 1, 00, 1, . . .), An is the set of all partitions of n into distinct parts, and 
Bn is the set of partitions of n into odd parts. The bijection i ^ 2i between supp(a) = P 
and supp(6) = 2P satisfies iai = so a b and \An\ = \Bn\- Throughout the 

paper, we refer to this example as the distinct/ odd case. 



2.2. O'Hara's algorithm. The analytic proof of Andrews's theorem shown above does 
not give an explicit bijection An — > Bn- Such a bijection is, by Theorem \'2.'6\ given by the 
following algorithm. 

Algorithm 2.2 (O'Hara's algorithm on partitions) 

Fix : sequences a ~^ b 
Input : X £ A 
Set : fj, <^ X 

While : /u contains more than bj copies of j for some j 

Do : remove bj copies of j from fi, add copies of i to fi, where (p{i) = j 
Output : ^{X) <— fi 

Theorem 2.3 (O'Hara) Alaorithm \2.2\ stops after a finite number of steps. The resulting 
partition ipiX) £ B is independent of the order of the parts removed and defines a size- 
preserving bijection A ^ B. 

Later on (see Subsection 14. 2p we deduce O'Hara's theorem from our generalization (Theo- 
rem [3]2]). 

Denote by L^{X) the number of steps O'Hara's algorithm takes to compute ipiX), and by 
Cip{n) the maximum value of L^{X) over all A K n. 

Example 2.4 In the distinct/odd case, O'Hara's algorithm gives the inverse of Glaisher's 
bijection, which maps A = I'^iS™'^ • • • £ B to the partition fj, £ A which contains if and 
only if rui has a 1 in the j-th position when written in binary. □ 
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Example 2.5 Let a = (1, 1, 4, 5, 3, 1, 1, . . .), ^ = (1, 1, 5, 3, 4, 1, 1, . . .) and ip{3) = 4, ip{4) = 
5, If (5) = 3, (p{i) = i for i / 3, 4, 5; observe that a b. Then O'Hara's algorithm on 
A = 3'^4^5^ runs as fohows: 

334452 _^ 374I52 _^ 324I55 ^ 324651 ^ 364351 
3104051 ^ 354054 _^ 304057 ^ 304553 ^ 3^425^ 

We have L<^(A) = £^(35) =9. □ 

Example 2.6 Take a = (2,2,1,2,2,1,...) and b = (3,1,3,1,...). Here A is the set of 
partitions into distinct parts = ±1 mod 3, and B is the set of partitions into odd parts, 
none appearing more than twice. Define ip: P — > P as follows: 

{i if i is divisible by 6 
i/3 if i is divisible by 3, but not by 2 . 
2i if i is not divisible by 3 

Clearly, a b. O'Hara's algorithm on 1^2^8^10^14^20^ runs as follows: 

1^2^8^10^14^20^ ii2i8iio3i4i 1I2I728I103 1I2I52728I102 

11215*728^10^ Ii2i5'^728i 1121425*^72 Ii234i5*'72 

II255672 _^ 132-15672 1^235672 1^225672 

192I5672 ^ 1II5672 _^ 1II537215I 1II72152 

^ 183I72152 ^ 153272152 123372152 ^ 127^9^15^ 

The bijection -0 is similar in spirit to Glaisher's bijection: given A = i™-i2"^'^4'^'^5"^^ ■ ■ ■ £ A 
and j G P, the number of copies of part 2j — 1 in tp{X) is equal to the k-th digit in the 
ternary expansion of I, where k is the highest power of 3 dividing 2j — 1, 2j — 1 = 3^r, and 

z = E^2W- □ 



2.3. Equivalent sequences and graphs. Choose equivalent sequences o, b. Define a 
directed graph on supp(a) Usupp(5) by drawing an edge from i to j if (p{j) = i; an arrow 
from i to j therefore means that O'Hara's algorithm simultaneously removes copies of i and 
adds copies of j. Each vertex v has indegv < 1, outdeg?; < 1 and indegf + outdegu > 1. 
The graph splits into connected components of the following five types: 

(i) cycles of length m > 1; 

(ii) paths of length m > 2; 

(iii) infinite paths with a starting point, but without an ending point; 

(iv) infinite paths with an ending point, but without a starting point; 

(v) infinite paths without a starting point or an ending point. 

Example 2.7 Figure [T] shows portions of graphs for certain ip: 

(1) a = (1,1,4,5,3,1,1,...), b = (1,1,5,3,4,1,1,...), ip{3) = 4, 99(4) = 5, ip{5) = 3, 
ip{i) = i for i 7^ 3,4, 5; components of G,p are of type ([i|); 

(2) a = (cxD, 1, 2, 3, 00, 00, 00, . . .), b = (2, 3, 4, 00, 00, 00, 00, . . .), ip{2) = 1, ip{3) = 2, 
ip{4:) = 3; G^ is of type (jn]); 

(3) the distinct/odd case: a = (2, 2, . . .), b = (00, 1, 00, 1, . . .), ip{i) = 2i; components of 
G^ are of type (fni]) : 
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(4) the odd/distinct case: a = (cxo, 1, oo, 1, . . .), 6 = (2, 2, . . .), ip{i) = i/2\ components 
of G^p are of type (IrvTl : 

(5) a = (2, 2, 1, 2, 2, 1, . . .) and & = (3, 1, 3, 1, . . .), if given by (|2.ip : components of 
are of types ([ij) and □ 
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Figure 1. Examples of graphs G^. 



2.4. Scissor-congruence and Il-congruence. We say that convex polytopes A, B in M*" 
are congruent, write Ac^ B,\i B can be obtained from A by rotation and translation. For 
convex polytopes P,Q C M™", we say that they are scissor-congruent if P can be cut into 
finitely many polytopes which can be rearranged and assembled into Q, i.e. if P and Q are 
the disjoint union of congruent polytopes: P = U^L^Pj, Q = U^^-^Qi, Pi ^ Qi. 

Let TT be a linear functional on M™'. If Qi can be obtained from Pi by a translation by a 
vector in the hyperplane 7i = {x £ R™: 7r(x) = 0}, we say that P and Q are ir-congruent. 
If P and Q are vr-congruent for some linear functional vr, we say that they are U-congruent. 

If P can be cut into countably many polytopes which can be translated by a vector in the 
hyperplane Ti. = {x £ : 7r(x) = 0} and assembled into Q, we say that P and Q are 
approximately TT-congruent. We say that they are approximately U-congruent if they are 
approximately vr-congruent for some linear functional vr. If P and Q are approximately 
TT-congruent, there exist, for every e > 0, vr-congruent polytopes Pe Q P and Qe C Q, such 
that vol(P \ Pe) < e and \o\{Q \ Qe) < e. 

Finally, let R(ai, . . . , am) = [0, ai) x • • • x [0, am) be a box in M™, and let R{ai, . . . , am) = 
R(ai, . . . , am) n Z™ be the set of its integer points. 

Example 2.8 Let d = 2 and vr(x,y) = x + y. Euclid's algorithm on (a, 6) yields a vr- 
congruence between R(a, h) and R(5, a): if 6 = ria+si with < si < a, divide [0, a) x [0, na) 
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into ri squares with side a, and translate the square [0, a) x [ia, {i + l)a) by the vector 
{ia,—ia) to [ia,{i + l)a) x [0, a). Then write a = r2Si + S2 with < S2 < si, divide 
[0, a) X [ria, b) into r2 squares with side si, and translate the square [isi, {i + l)si) x [ria, h) 
by the vector (rio — isi^isi — rio) to [rio, 6) x [isi, (i + Continue until the remainder 

Si is equal to 0. The first drawing of Figure [2] gives an example. 

The second drawing shows that boxes R(12, 8) and R(32, 3) are vr-congruent for 7r(x, y) = 
X + 4y. Finally, in Figure [3] we give a vr-congruence between R(4, 5, 3) and R(5, 3, 4) for 
7r(x, y, = 3j; + 4y + Sz. □ 




1 4 




^ 1 3 








3 1 4' 


2' 



Figure 2. Two H-congruences. 





Figure 3. vr-congruence between R(4, 5,3) and R(5,3,4). 



3. Main results 

3.1. Continuous O'Hara's algorithm and H-congruences. Take the case when G^p 
is a cycle ii ^ im ^ im-i ii. In this case, 97(^1) = 12, ^{12) = ^3, etc. 



Throughout this section, identify a partition 



with the vector t 



By Theorem 12.31 O'Hara's algorithm defines a bijection ip: i?(ai, . . . ^cirn 
where ijaj = ij^ibj^i for all j. The following algorithm (see also Theorem 13. 2p general- 
izes ijj to the continuous setting. It gives a bijection i/> : R(ai, . . . ,am) — >■ ^{bi, . . . ,bm), 
which is defined also for non- integer aj,bj. When aj,bj are integers, it is an extension of 
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ip: R{ai, . . . ,am) ^ ■ • • , &m)- As an immediate corollary, we prove that two boxes with 
rational coordinates and with equal volume are Il-congruent. We show in Subsection 14.21 
how we can use Theorem 13.21 to give an alternative proof of Theorem 12.31 

Algorithm 3.1 (continuous O'Hara's algorithm) 

Fix: i = {h,...,im)eR^ 

a = (ai, . . . , am) G W^, h = {bi, . . . ,bm) ^ with ijaj = ij+ibj+i 
Input : t G R(ai , . . . , am) 
Set : s ^ t 

While : s contains a coordinate Sj > bj 

Do : Sj ^ Sj — bj, Sj_i <— Sj-i + aj-i 
Output : ■j/j(t) <— s 

It is clear that the algorithm starts with an element of P = R(ai, . . . , am) and, if the while 
loop terminates, outputs an element of Q = R(5i, . . . , bm)- It is not obvious, however, that 
the loop terminates in every case, or that the output i/'(t) and the number of steps L^(t) 
depend only on t, not on the choices made in the while loop. 

Theorem 3.2 Algorithm \3. 1\ has the following properties. 

(1) The algorithm stops after a finite number of steps, and the resulting vector xj^{t) and 
the number of steps L^(t) are independent of the choices made during the execution 
of the algorithm. 

(2) The algorithm defines a bisection if^: P ^ Q which satisfies i/'(t) — t G 7i, where 7i 
is the hyperplane defined by iixi + . . . + imXm = 0. 

(3) We have 

L^(t + t') > L^(t) +L^(t') for every t,t',t + t' G P. 

In particular, L^(t') < L(p(t) if t' < t. 

(4) Let t, t' G P, s = ^l){t), with tj < t'j < tj + Sj, where Sj = bj — Sj. Then 

t/j(t') - t' = t/>(t) - t and L<^(t') = L^(t). 

(5) For all a, b G Z™, we have 

maxL;p(t) = Icmfci, . . . ,Cm) ■ ( — + • ■ • + — ) - rn, 

tGP ^ \Ci Cm) 

where Cj = a\ - • • aj-\bj ■ ■ ■ bm-i ■ 
The proof of the theorem is given in Subsection 14.11 

We call boxes P = R(ai, . . . , Om), Q = R(^i, • • • , bm) relatively rational if there exists A, 
A 7^ 0, such that Xoj G Xbj G Z. Clearly, two boxes P and Q with rational side-lengths 
are relatively rational. 

Corollary 3.3 Boxes P = R(ai, . . . , a^), Q = R(6i,...,6m) with equal volume are ap- 
proximately li-congruent. Moreover, when P and Q are relatively rational and have equal 
volume, they are U-congruent. 
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Proof. For j = 1, . . . , m, take ij = ai • • • aj_i5j+i • • • bm- Clearly ijaj = ij+ibj+i for j = 
1, . . . , m — 1, and ai ■ ■ ■ am = bi • • • bm implies imam = ^i^i- Therefore, the numbers ij, aj, bj 
satisfy the conditions of Algorithm 13.11 By Theorem 13.21 part ([2]) , the algorithm defines 
a bijection ■0: P — > Q. Parts ([4]) and ([2]) of Theorem 13.21 imply that we can cut P into 
(countably many) smaller boxes, each of which is translated by a vector in the plane iixi + 

If P and Q are relatively rational, we can assume without loss of generality that all aj , bj 
are integers. For any integer vector t, we have ■j/'(t') — t' = ■j/'(t) — t and L^(t') = L^(t) 
whenever tj < t'- < + 1, so P and Q are divided into a finite number (at most oi • • • Om) 
of boxes. □ 

Example 3.4 Even in the 3-dimensional case the H-congruence defined by the algorithm 
can be quite complex, as the next figure suggests. Here the same shading is used for parallel 
translations by the same vector. □ 




Figure 4. The decomposition of the box R(31,47, 23) given by O'Hara's 
algorithm (only the top, right, and back sides are shown) . 



3.2. Complexity of O'Hara's algorithm. The complexity of O'Hara's algorithm has 
been an open problem, with the exception of the elementary distinct /odd case (see [01] 
and Example 13. Sh . 

It turns out that the complexity depends heavily on the type of the graph dp defined in 
Subsection 12.31 Part ([5]) of Theorem 13.21 gives the maximum number of steps that O'Hara's 
algorithm takes when is a cycle. The following lemma gives an estimate for Cp^n) when 
Gp is a path. 
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Lemma 3.5 Let be a finite or infinite path on I Q ¥. Then Cip{n) < n(logn + 1). 
Moreover, if 



then C^{n) < Dn. 



This lemma and the other results in this subsection are proved in Section [4l Here and 
throughout the paper, by logn we mean the natural logarithm of n. 

We combine these estimates to prove Theorem 13.61 the second main result of this paper, 
which gives upper bounds for the complexity of the algorithm in the general case. 

Theorem 3.6 Let a, b be ip-equivalent sequences. 

(1) If has only a finite number of cycles of length > 2, then Cip{n) = 0{nlogn), 
and the constants implied by the O-notation are universal. 

(2) // Gip has only a finite number of cycles of length > m for some m > 2, then 
■-ipi 



C^p{n) = 0{n"^ and the constants implied by the O-notation depend only on m. 



The following theorem gives the corresponding lower bound on the worst case complexity. 
It shows that the estimates of Theorem 13.61 are close to being sharp. 

Theorem 3.7 There exist ip-equivalent sequences a and b, such that: 

(1) G^p is a path and C^p{n) = log logn); 

(2) G^p contains only cycles of length < m and C^{n) = n(n™'~^~^) for every e > 0; 

(3) £^(n) = expJ7(^). 

In other words, depending on the type of the graph, we have nearly matching upper and 
lower bounds on Cip{n). For example, for an m-cycle, Theorem 13.61 shows that Cip{n) 
is 0{n^~^), while Theorem 13.71 shows that it is Q{n^~^~'^) for every e > 0. Similarly, 
part ([3]) shows that O'Hara's algorithm can be very slow in general since the total number 
of partitions of n is asymptotically exp0(-y/n). 

Example 3.8 In the distinct/odd case, the graph G^ is composed of infinite paths 

. . . — > 8j ^ 4j ^ 2j j for each odd j. 

A partition X h n can be broken up into partitions A'^-'^ h uj such that the support of A*^-'^ 
is contained in {(2j — 1)2^^ : k G N}. We have 



oo ^ 



^^^2^+i(2j-l) 2j-l' 

and Lemma [3.51 implies that O'Hara's algorithm takes at most nj/{2j — 1) steps to compute 
^'(A^-'^). This implies that 

oo n 
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In other words, O'Hara's algorithm takes at most n steps to compute ip{X)- This bound is 
(almost) sharp since the algorithm takes 2^ — 1 steps to compute 11^(2^ ) = 1^ . □ 

3.3. O'Hara's algorithm as an integer linear programming problem. Let us now 

give a new description of O'Hara's algorithm. 

Proposition 3.9 Let i, a, b G be as above such that ijUj = ij+iftj+i for j = 1, . . . , m. Fix 
a vector t G R(ai, . . . , am)- Then s = i/'(t) satisfies the following: 

s = t + Ak, 

where 



f-b. 


ai 


•• 


• ^ 





-b2 


a2 ■ ■ 











-bs ■■ 










•• 


• -bmj 



and k = {ki, . . . , km) is the unique vector minimizing 

^1 + . • • ~l" kfii 

with constraints 

k G Z"", k > 0, ^k > -t, ^k < b - 1 - t. 

The proposition is proved in Subsection 14.61 The advantage of this approach is that one 
can use standard integer linear programming results to speed up the computation of i]j{t). 

Theorem 3.10 For every m > 1, there exists a deterministic algorithm which computes 
the continuous O'Hara's bijection ip in polynomial time, for all integer vectors i,a, b as 
above. 

Proof It is well known that for a bounded dimension m, there exists an algorithm for solving 
integer linear programming problem Ax < b, for which the number of steps is bounded by 
a polynomial in the logarithm of the largest entry of A, h for integer A and b (see e.g. [HI 
Corollary 18.7b]). By Proposition 13.91 this implies the result. □ 

Theorem 13 . 1 1 can be used to obtain a significant speed-up of (the usual) O'Hara's algorithm, 
in the case when contains only cycles of bounded length. Namely, we obtain the following 
result. 

Theorem 3.11 Let a ~^ b. If the lengths of cycles of are bounded, there exists a 
deterministic algorithm which computes ip{X) in O(nlogn) steps for X G An- 

Proof. Without loss of generality, the support of A £ An is contained in one of the con- 
nected components of G^. If this connected component is a path, O'Hara's algorithm takes 
0(n log n) steps by Lemma 13.51 If it is a cycle of length m, the algorithm described in 
the previous theorem takes O(log'^n) steps for some c, and obviously the 0(n log n) term 
dominates. □ 
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Remark 3.12 Let us note that the inner workings of the algorithms in Theorem 13.101 and 
Theorem 13.111 have a geometric rather than combinatorial nature, and are very different 
from those of O'Hara's algorithm. However, both kinds of algorithms, when applied to 
the same input, have the same output, which means that they produce the same partition 
bijection. 



4. Proofs and examples 
4.1. Proof of Theorem 13.21 Throughout the section, indices are taken modulo m. 

Lemma 4.1 Take a vector t £ P, choose £j < aj — tj, and do the algorithm on t: denote 
the vectors we get by = t, , s^, etc. Then each box of size ei x ■ ■ ■ x Em contains at 
most one of s* . 

Proof. Assume that we have |s^' — s^l < ej for all j for some < k < k' , i.e. that we hit 
an El X ■■■ X Em box twice. Say that in the course of getting from t to s*^ (respectively 
s'^ , respectively s'^ we subtracted bj from the j-th. coordinate and added aj to the 
(j — l)-th coordinate kj times (respectively k^ times, respectively kj times). Clearly we have 
kj = k,Y^ kj = k' and ^ k'- = k' — k. Furthermore, the equations Sj = tj — bjkj + ajkj+i, 

Sj' = tj — bjk'j + ajkj_^^ and = tj — bjk'- + ajkj_^i, which hold for all j, can be written 
as 

s*^ = t + Ak, 
8^=' = t+^k', 
s'='-'= = t + Ak", 

where 





f-bi 


ai 


•• 


• ^ 







-62 


02 • • 





A = 








-63 •• 












• • 


• -bmj 



There are two cases to consider: k" = k' — k and k" 7^ k' — k. We obtain a contradiction 
in each case. Assume first that k" = k' — k. Then we have 

s^''-'^ - t = t + A{k' - k) - t = s'^' - 

In particular, we have \sj ~^ — tj\ < Ej and Sj ~^ < Oj for all j. On the other hand, s'^'"'^ was 
obtained from s'^'~'^~^ by choosing j with Sj "'^"^ > bj and then taking s^ ~^ = -'^-i 
and s'^'Si'' = s'^'Si^^ + Oj-i- Therefore Sj'j^'' > Oj-i, a contradiction. 

Suppose now that k" 7^ k' — k. Since ^j — Yli^ ~ there is a j so that kj > kj — kj, 
k'j+i < ^j+i ~ ^j+i- We have 

s'^'-'^ = t + Ak!' = t + A{k!' - k' + k) + 8^=' - 
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This implies 

< s^-~^ < tj - bj {k'l - k'j + kj) +aj (fcj+i - A;^+i + kj+i) +ej < tj - aj + £j < 0, 

^ V ' ^ V ' 

>0 <-l 

a contradiction. Therefore we can never hit an ei x • • • x box twice, which completes 
the proof. □ 

Now we are ready to prove all parts of Theorem 13. 2i 

JH) Clearly we have iis^ + . . . + ims'^ = iih + . . . + imtm and Sj>0 for all j and k. Since 
the set {x G : iixi + . . . + imXm = C, xj > 0} can be covered by a finite number of boxes 
of size ei X • • • X e^, the number of steps is finite. 

Let us prove that the output vector s £ Q does not depend on the choices we make in the 
course of O'Hara's algorithm. Assume that when we run the algorithm twice, we obtain 
vectors s and s'. Denote hy k = |k| = ki + . . . + km the number of steps to obtain s, and 
by k' = |k'| the number of steps to obtain s', and assume that k < k' . Denote by s" the 
vector we obtain after k steps in the second run of the algorithm. Think of s" as being 
on the path from t to s'. Write s = t + Ak, s" = t + Ak", where |k| = |k"| = k. Then 
s" = s + A(k" - k), and kj > kj, kj^^ < kj^i would imply 

< s'j = Sj — bj {kj — kj) +aj {k'-^i — /cj+i) < bj — bj = 0, 

' >i ' ' <o ' 
which is a contradiction. Therefore k" = k and s" = s, therefore also k' = k and s' = s. 

(I2|) Let us construct an explicit inverse map jjj ^. Denote by t the flip t: (xi, . . . , Xm) ' — ^ 
{xm, ■ ■ ■ ,xi). Consider a map : R(6m, • • • , ^i) — > R-lom, •••,01) obtained by replacing 

[ll , ■ ■ • , 1m 

), (ai,...,am) and {bi,...,bm) with {im,---,ii), {bm,---,bi) and ( 
respectively. Observe that 

^-1 :=ro'0'or:Q P 
is the inverse map of tp. Therefore, ■0 is one-to-one. 

It remains to check that ^p satisfies ip{t) — t £ TC, where Ti. is the hyperplane defined by 
iixi + . . . + imXm = 0. Note that the columns of A lie in TC. Because i/'(t) — t = Ak for 
some k, ip{t) — t lies in TC. 

([3]) Let = t,s-'^,s^, ... ,8*^ = tp{t), k = L(^(t), be the intermediate steps of the algorithm 
which computes ■0(t). Similarly, let r° = t', r-*^, r^, . . . , r*^' = i/'(t'), k' = L(^(t'), be the 
intermediate steps of the algorithm which computes i/'(t). Every vector in the sequence 
t + t', + t', . • • ) s'^ + t', s*^ + r^, . . . , s'^ + r'^' ~^ has at least one of the coordinates > bj, so 
O'Hara's algorithm takes at least k + k' = L^(t) + L^(t') steps to compute ■?/'(t + 1'). 

(jl]) Let t G P, s = ■j/'(t) and Sj = bj — Sj. For every t' with tj < t'j < tj + £j, make the steps 
of the algorithm which inputs t' the same as the one which inputs t. Let s* and r* be as in 
part dH]). Then the vector r* satisfies r* = s* + (t' — t). For k = L^(t), we have s'^ = s and 
= s + (f - t) G Q. Therefore, r'^' = 1/7(1'), ^(f) = i/7(t) + t' - t and L<^(t') = L^(t), as 
desired. 
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(j5]) Assume that aj and bj are integers for all j = 1, . . . ,m. Let t £ P and assume that 
O'Hara's algorithm with input t subtracts bj from the j-th coordinate and adds flj-i to the 
{j — l)-st coordinate exactly kj times. This implies that s = i/'(t) = t + ^k. 

Observe that the matrix A has rank m — 1, and its kernel is spanned by the vector 
1 1 1 \ / 1 1 1 



^6l52 • • • &m-2^m-l ^1^2 ' ' ' ^m-2&m-l «1«2 ' ' ' am-2am-l / V^l C2 C„ 

The integer vectors in ker A are integer multiples of the vector 

lcm(ci,C2,...,Cm) • ( — , 

Suppose an integer vector n = (ni,...,nm) with at least one negative and at least one 
non-negative coordinate is a solution of 

(4.1) s = t + An. 

Suppose rij < and ?t-j+i > 0. Then we have —bjUj + ajUj+i > bj and 

bj > Sj = tj — bjUj + QjUj-^i > bj, 

a contradiction. This implies that the coordinates of every solution to (j4.ip are either all 
negative or all non-negative. In particular, there is exactly one solution n with 

< n, < - 1 

for each j. 

First, let us show that this solution is equal to k. Denote by n'j the number of times 
we subtract bj from the j-th coordinate and add Oj-i to the {j — l)-th coordinate in the 
first |n| = 77-1 -|- ... -|- rim steps of the algorithm with input t. Let s' = t -|- An'. Then 
s' = s -|- ^(n' — n), and the same argument as above shows that we cannot have n'j — rij < 0, 
n'j_^_'^ — rijj^i > 0. Clearly the sum of the coordinates of n' — n is 0, so we must have n' = n 
and s' = s. This implies that n = k. Therefore, the maximum number of steps L^(t) is at 
most 

(4.2) lcm(ci,...,Cm) • {— + ... + —) -m. 

Consider the vector (ai — 1, . . . , — 1) and observe that 

n(ai - 1) -h . . . imiam - 1) = n(&i - 1) + . . . + im(&m " !)• 

Since i/'(t) is an integer vector, we have s = i/'(t) = (6i — 1, . . . , 6^ — 1). Furthermore, the 
vector (— 1,— 1,...,— l)isa solution of the system 

bj — \ = aj — \ — bjUj + ttjUj^i, 1 < j < Tn. 



Therefore, 



_ lcm(ci, . . . ,Cn^) 



and the maximum number of steps L(^(t) over all t, is given by the equation (14. 2p . This 
finishes the proof of Theorem 13. 2[ □ 
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4.2. Proof of Theorem l2.3l In this subsection, we present an alternative proof of O'Hara's 
theorem by deducing it from Theorem 13.21 Take a partition A G An- Without loss of gener- 
ality we may assume that its support is contained in one of the connected components of G^. 
If this connected component is a cycle, the fact that O'Hara's algorithm stops and that the 
final result is independent of the choices made in the algorithm follows from Theorem 13.21 
part ([T]). If the connected component is a path, we can assume that A = • • • i^, where 
ii ^ 12 ^ ■ ■ ■ ^ im ^ im+1 ^ ... is a part of G^p. Assume that ijhj = ij+iOj+i. Remove 
copies of ii until their number is smaller than bi. Remove copies of 12 until their number 
is smaller than b2, etc. Since ij are all different, we have ik > n for some k. The algorithm 
cannot increase the number of copies of it, so it stops before that. Furthermore, it is clear 
that if kj is the number of copies of ij removed, then ki = k2 = [(t2 + fcia2)/^2j7 

etc. Therefore the final result is independent of the order in which we remove and add 
parts. □ 



4.3. Proof of Lemma 13. 5L Assume that the path is ... i_i — > zq — > ii — > 22 — *■ • • • 

(it can be finite, or infinite in one or both direction). A partition A has only a finite 
number of parts, without loss of generality we can assume that A = i*ii^2 ' ' ' ^ni '^{^) — 
i^^ ■ ■ ■ if^ . O'Hara's algorithm is straightforward in this case: first remove biki copies 
of ii, where ki is the smallest integer with ti — biki < bi, and add 02/^1 copies of 12; then 
remove 62/^2 copies of 12, where /c2 is the smallest integer with t2 + 02/^1 — 62^2 < &2) and 
add 03/02 copies of is, etc. Clearly, 

, . ^1 kh 
, ^t2 + a2ki i2t2 + kbiki ^ kh + kh 

^2 < 7 = TT < TT , 

O2 Z202 ^202 

, + 03^2 ^3*3 + i2&2fc2 , ^1*1 + ^2*2 + ^3*3 

^3 ^ 1 = -I ^ ' 

O3 ^303 «303 

etc. Therefore, the total number of steps is at most 

{kti + ...+ imtm) ( + . . . + 



kbi imbn 

Since kk + . . . + imtm is the size n of A, the number of steps is at most Dn when D < oo. 
In general, since ii, . . . ,im are distinct integers, there is a j with ij > m. On the other 
hand, when ij > n, O'Hara's algorithm stops at ij. This implies that m <n, and 

1 11 111, 

+ . . . + ^-r- <— + ... + — <1 + - + ...H <logm + l<logn+l. 



iih imbm ii im 2 m 

This completes the proof. □ 
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4.4. Proof of Theorem 13. 6L (1) For a 1-cycle {i} and a partition A with only parts i, 
we have ip{X) = A and O'Hara's algorithm takes steps. For a 2-cycle i — > j — > z, by 
Theorem 13.21 part ([5]), the largest number of steps of O'Hara's algorithm is equal to 

lcm{ai,bi)- ( — + r ) - 2<ai6i - ( — + ^ ) - 2 = + 6^ - 2. 
V«i hj \ai bij 

Denote by Amax = the largest partition in A. Then O'Hara's algorithm takes at 

most 

ai + bi-2 < i{ai - 1) + i{bi - 1) = |Amax| steps. 

For a partition A h n > |Amax|/2, we need, by Theorem 13.21 part ([3]), at most |Amax| < 2n 
steps. Similarly, for a partition A = h n < |Amax|/2, we have L^{X) < max{t,s} < n, 
since it + js is less than — 1) or j{bj — 1). 

Consider a cycle C : n ^ im — > im~i ^ . . . — > ii in dp of length m > 2. Denote by Mc 
the maximum number of steps O'Hara's algorithm can take on partitions with support in 
{ii, . . . ,im}- Recall that Mc is given by Theorem 13.21 part ([5]). Finally, denote by M the 
sum of Mc over all cycles C of length > 2 in G^. 

Note that every partition A E ^ is decomposed into partitions \^^\ . . . , X^"^^ of sizes ni, . . . , 
with support in one of the components of G^. From above and by Lemma 13.51 the number 
of steps of O'Hara's algorithm is at most 

ni(logni + 1) + . . . + nr{lognr + 1) + 2n + M < n(logn + 3) + M < (1 + e)nlogn, 

for every fixed e > and n large enough. 

(2) Take a cycle ii — > — *■ im-i of length m. In the notation of Theorem 13.21 

part ([5]), we have 

L^{X) < lcm(ci, . . . ,Cm) ■ (7- + . • • + — ) - m 



< 0102 • ■■ara-\b\b2 ■ ■ ■ bm-1 ■ ( — + • • • + — ) - m 

\Cl Cm) 

m 

- ^1^2 • • • bj-iaj ■ ■ ■ am~2am-i - m. 
i=i 

Denote by N is the largest of {oi, . . . , Om, 61, ... , bm}- Then the size of the largest partition 
Amax G .4. is at least — 1, and by the calculation above the number of steps is at most 
m{N"^~^ — 1). Therefore L^(Amax) < Amaxl™""^- For A h n > |Amax|/"T.i we have 
Lip{X) < ?Ti| Amaxl™'"^ ^ m™'n'"~^. On the other hand, when n < \Xaia.x\/'fn, we have 
n < ij{bj — 1) for some j, and O'Hara's algorithm does not remove any copies of ij. By 
Lemma 13.51 this implies that Lip{X) < n{l + 1/2 + . . . + 1/m). 

In the notation of part (1), write M' for the sum of Mc over all cycles of length > m. As 
before, every partition A € An is decomposed into partitions A'^^^ . . . , A^*") of sizes ni, . . . , 
with support in one of the components of G^. From above and by Lemma 13.51 t^is number 
of steps of O'Hara's algorithm is at most 

^m^m-i + _ , + m™n;"-^ + M' < m^n*"-^ + M' = 0(n"-i). □ 
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4.5. Proof of Theorem 13.71 We begin with the following key example. 

Example 4.2 Let pi,P2, ■ ■ ■ ,Pm be distinct primes. Take ip{pj) = Pj+i (where Pm+i = Pi), 
aj = pj+i, bj = pj-i (where po = Pm)- Then 

0102 • • • aj_i6j • • • 6m_2&m-l = P2P-i ' ' ' Pj-lPjPj~lPj ' ' ' Pm-2 

for j = 1, . . . ,m, and the lowest common multiple of these numbers is 

PlPhl ■ ■ ■ Pm-2Pm-lPm- 

By Theorem 13.21 part ([5j), the maximum number of steps of O'Hara's algorithm is equal to 

/II 1 1 \ 
PiP2 ■■■Pm- \ h . . . H \ m. 

\PlP2 P2P3 Pm-lPm PmPl J 

Note that Example 12.51 is a special case of this for m = 3, pi = 3, p2 = 4 and p^ = 5. 

When are approximately equal, i.e. \pj/N — 1| < e for some N, we obtain 

partitions of relatively small size (around mN"^), for which O'Hara's algorithm takes a large 
number of steps (around mN^~'^). □ 

Now we are ready to prove Theorem 13.71 
(1) Take 



^3 


1 - 


2 - 


3 - 


6 - 


5 - 


10 - 


7 - 


14 - 


9 ^ ... 




oo 


1 


2 


1 


6 


1 


10 


1 


14 




2 


3 


2 


5 


2 


7 


2 


9 


2 



i.e. i2j-i = 2j - 1, i2j = 4j - 2, ip{ij) = ij-i, ai = oo, 02^-1 = 4 j - 6 for j > 2, 04^-2 = 1, 
62j-i = 2, 64^-2 = 2j + 1. O'Hara's algorithm on the partition A'-'^^ = 1^ (sfc-i)!! j-^ng as 
follows: 



-j^2''-3-5-(2k-l) 

22*~l-3-5---(2fc-l) _ 


22*^ 


"l-3-5---{2fc-l) 
l-5-7----(2fe-l) 


in 


2k- 


1 


3 


5 


7 


9- 


■{2k- 


1) 


steps 


32'=-l-5-7----(2fe-l) _ 
g2'=~l.l-5-7---(2fe-l) _ 


62*" 
-y 52'= 


^■l-5-7---(2fc-l) 
■l-3-7---(2fc-l) 


in 


2k- 


1 


1 


5 


7 


9- 


■ {2k- 


1) 


steps 


52'=-l-3-7---(2fc-l) _ 


102'=" 


-1 .1.3.7. ..(2fc_l) 


in 


2k- 


1 


1 


3 


7 


9- 


■{2k- 


1) 


steps 



Note that n := |A(^)| = 2*=(2/fc - 1)!! = k'^^'^l We conclude that 

> 2' . (2k - 1)H (1 + 1 + 1 + , , , + ^pi-^) = !i(„l„gl„g„). 

(2) Take to be a superpolynomial (for example exponential) integer function of k. For 
each k large enough, we can choose m distinct primes Pi, ■ ■ ■ ,Pm between Dj^ and 20^. 
Choose ij = km + j for j = 1, . . . , m, = ij^iPj, bj = ij_iPj_i (where indices are written 
cyclically). Therefore, ij^ — > ^ . . . ^ is an m-cycle in G^. The largest partition 
A^*') G A with support in {i^, . . . , i^} has size at most 2{k + \)m?Dk- Then 

lcm(ci, . . . , Cm) = \cm{PmPlP2 ■ ■■Pm-2,P\P2 ' ' ' Pm-2, ■ ■ ■ ,PlP2 ' ' ' Pm-l) 
= P1P2 ■ ■■Pm-2Pm-lPm 
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where we omitted the upper index k of for simphcity. Therefore, 

L^{\^^'^) > V1P2 ■ ■■Pm-l +P2P3 ■■■Pm + ■■ ■+PmPl ' ' ' Pm~2 -m> m{D'^~^ - 1). 

On the other hand, for every C and k large enough, we have 

m{D^-^ - 1) > C{2{k + iWDkT'^-', 

since otherwise grows polynomially, contradicting our assumptions. This imphes that 
the corresponding a, b and ip satisfy the conditions of (2). 

(3) Let pi,P2, ■ ■ ■ denote the sequence of ah primes. Set (p{pj) = Pj+i for j / /c'^, ip{pi) = pi, 
= P(k~i)''-^+i for k>2, aj = pj+i for j / k'', ai = pi, a^k = for k > 2, 

bj = Pj-i for j / 1, j / A;'^ + 1, bi = pi, ^^fc+i = P(fc+i)fe+i for k > 1. The following table 
summarizes these values. 
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7 


11 


13 


17 
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23 
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31 


37 


41 


43 
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103 
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47 
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2 


7 


3 


5 


103 


11 


13 


17 


19 


23 


29 


31 


37 


41 


43 . 


. 97 


101 





Clearly, dp consists of cycles of length 1 and k^ — {k — 1)^ ^ for k > 2. For k > 2, write 
mi = {k — l)'^"^, m2 = k^ , and define A^'^^ S ^ to be the partition 



of size 

nk = Pmi+lPmi+2 + Pmi+2Pmi+3 + • • • + Pm2-lPm2 + Pm2Pmi+l — Pmi+l — Pmi+2 — ■ ■ ■ — Prn.2 • 

By the calculation in Example 14.21 O'Hara's algorithm takes exactly 

Pmi+lPmi+2 ■ ■ ■ Pm2-2 + Pmi+2Pmi+3 ' ' ' Pma-l + • • • + PmaPmi+l " " " Pma-S " " "T-i) 

steps to compute By the distribution law of prime numbers, we have pn = 

nlogn(l + 0(1)). Therefore, 

nk < (m2 - m^)pl^ = {k^ - {k - lf-^)pl, ~ log(A;'=))2 = k''''^^ {\ogkf 

and 

logL^(AW) > log ((m2 - mi)(p^^;fi-2 - 1)) ~ kHog {{k - if-^) ~ A^'^+MogA:. 
Thus, LpiX^^'^) > exp(^/7^) for k large enough, as desired. □ 

4.6. Proof of Proposition [3T9l By construction of the algorithm, k G and it satisfies 
the inequalities in Proposition 13.91 Therefore, it suffices to prove that |k| < |n| for every 
non-negative integer vector n 7^ k satisfying < t + An < b — 1 . 

Assume that |n| < |k|. Denote by n'j the number of times we remove bj copies of ij in the 
first |n| steps of O'Hara's algorithm. Define s' = t + An' . From above, s' = s + A{n' — n). 
If there exists j such that n'j — nj > and "-j+i — ^j+i < 0, then 

< s^- = Sj - bj{n'j - TLj) + aj(n^+^ - nj+i) < bj + {-bj) = 0, 
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which is impossible. Therefore, all the coordinates of n' — n have the same sign. On the 
other hand, the sum of the coordinates of n' — n is 0, which implies n' = n = k. Therefore, 
for all n 7^ k as above, we have |n| > |k|, as desired. □ 



5. Application: the speedy O'Hara's algorithm 

In \P3\ §8.2], a simple speed-up of O'Hara's algorithm was given: at each step, if the number 
of parts i in A h n is rbi + s for < s < 6j, remove rbi copies of part i and add raj copies 
of part j, where = i- This replaces r steps of O'Hara's algorithm with one step. This 
algorithm produces the same output as the original algorithm, but the number of steps 
may depend on the choices we make in the execution. It is called the speedy O'Hara's 
algorithm. Observe now that it cannot be much faster than the usual O'Hara's algorithm: 
since r < i{rbi + s) < n, the number of steps can be reduced by at most the order of 
n, where X h n. Thus the reasoning from the proof of Theorem 13.71 part ([3]) still gives 
superpolynomial lower bounds for the number of steps. The following examples show that 
the speed-up is not substantial (for example, of logarithmic complexity) even in the case 
described in Theorem l3.6l part ([T]), and that it gives the same bound as Theorem 13 . 61 part ([2]). 

Example 5.1 Take 







^ 18 




14 ^ 


5 


^ 10 - 


3 ^ 


6 


1 






1 


18 


1 


14 


1 


10 


1 


6 






7 


2 


5 


2 


3 


2 


1 


oo 


-1 = 


2j- 


1, hi = 


4j + 2, (^(ij) 


= ij+i, 


a2j_i 


= 4j + 2, 


a4j+2 = 1, 




= oo, 62 j 



i.e. = 2j - 1, i2j = 4j + 2, ip{ij) = ij+i, a2j-i = 4j + 2, 04^+2 = 1, 61 = 00, 62^-1 = 2 
for j > 2, 64j_|_2 = 2j — 1. This is similar to the example in the proof of Theorem 13.71 
part ([1]). The speedy O'Hara's algorithm on A*^'^) = {2k — l)4fc-2 £qj. ^ > 5 runs as follows: 

(2k - l)4k-2 ^ (4fc_2)2fc-i ^ (4fc-2)2(2fc-3)*'=-2 _ (4fc-2)2(4fc-6)2'=-i ~* (4fc-2)2(4fc-6)4(2fc-5)*'=-« 

^ (4A;-2)2(4fc-6)-4(4A:- 10)2^-3 ^ (4A:-2)2(4fc-6)4(4A:-10)4(2fc-7)-4'=-i0 ... ^ (4fc-2)2(4A:-6)'*--14'*10'^ 
^ (4fc-2)2(4fc-6)'*-144l0i320 -» (4fc-2)2(4fc-6)*-14*10i6io ^ (4k - 2)^(4k - 6)^ • • • 14^10^1*'°. 

In particular, the speedy O'Hara's algorithm takes 2k — 2 = Q{^/n) steps for a partition of 
size n = (2/c - 1)(4A: - 2) = 2(2A; - 1)2. □ 

Note that the proof that m < n in the proof of Lemma 13.51 together with the proof of 
Theorem 13.61 part ([T|) shows that the speedy O'Hara algorithm takes 0{n) steps when 
has a finite number of cycles of length greater than 2. This is smaller than r2(n log log n) 
obtained earlier. In other words, the speedy O'Hara's algorithm can be asymptotically 
faster. 

Example 5.2 Take a, b and constructed in the proof of part ([2]) of Theorem 13.71 see 
Subsection 14.51 For each k, we have i^/^i < {k + l)/k < 2, so at each step the number of 
parts ij is at most 

1 (iliilp'l - 1) + . . . + zH^pI - 1)) < Smi'lpl 
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Therefore, in the speedy algorithm we replace bj = ij_iPj_i copies of ij by Uj copies of 
ij-i at most 8m times. This means a speed-up by only a constant factor, and the number 
of steps of the speedy O'Hara's algorithm is 17(72™"^"^) for every e > 0. □ 



6. Final remarks 

6.1. The polynomial time algorithm in the proof of Theorem 13.101 is given implicitly, by 
using the general results in integer linear programming. It is saying that the function 

: An — > Bn can be computed much faster, by circumventing the elegant construction 
of O'Hara's algorithm. It would be interesting to give an explicit construction of such an 
algorithm. 

In a different direction, it might prove useful to restate other involution principle bijections 
in the language of linear programming, such as the Roger s-Ramanuj an bijection in [GM2] 
or in jBPj . If this works, this might lead to a new type of a bijection between these two 
classes of partitions. Alternatively, this might resolve the conjecture by the second author 
on the mildly exponential complexity of Garsia-Milne's Rogers- Ramanuj an bijection, see 
[P3l Conjecture 8.5]. 

6.2. Note the gap between the number exp0(-^/n) of partitions of n and the lower bound 
Cip{n) = expO(-y/n) in Theorem 13.71 It would be interesting to decide which of the two 
worst complexity bounds on the number of steps of O'Hara's algorithm is closer to the 
truth. 

Note that we applied our linear programming approach only in the bounded cycle case. We 
do not know if there is a way to apply the same technique to the general case. However, we 
believe that there are number theoretic obstacles preventing that and in fact, computing 
O'Hara's bijection as a function on partitions may be hard in the formal complexity sense. 

6.3. It would be interesting to find ip such that the graph is a path and C^{n) = 
0(77. log n). From the proof of Lemma 13.51 it follows that the number of steps of (the usual) 
O'Hara's algorithm is at most n(logA; + 1), where k is the number of steps of the speedy 
O'Hara's algorithm. In Subsection [5] we constructed an example with k ~ so it 
certainly seems possible that such examples exist. 

6.4. Most recently, variations on the O'Hara's bijection and applications of rewrite systems 
were found in |SSMj and \K1\ IK2j . It would be interesting to see connections between our 
analysis and this work. 

6.5. In the finite dimensional case, the structure of the map ■0 establishing approximate 
H-congruence remains largely unexplored. For example, it would be nice to obtain some 
convergence result in the irrational case using the rational approximations which follow 
from part (3) of Theorem 13.21 and upper bounds in part (2) of Theorem 13.61 

Recall also that the 2-dimensional case can be viewed as the Euclid algorithm which in turn 
corresponds to the usual continued fractions (see Example 12. 8p . Thus the geometry of t/' 
can be viewed as a delicate multidimensional extension of continued fractions. Given the 
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wide variety of (different) multidimensional continued fractions available in the literature, 
it would be interesting to see if there is a connection to at least one of these notions. 
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